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WIND-TUNNEL INVESTIGATION AT LOW SPEED OF THE AERODYNAMIC 
CHARACTERISTICS OF THE AH'^Y CHEMICAL CORPS MODEL 
E112 BOMBLETS 

By William Letko and James L. Williams 


SUMMARY 


An experimental investigation has been made in the Langley stability 
tunnel to determine the aerodynamic characteristics of the Array Chemical 
Corps model E112 bomblets. A detailed analysis has not been made; how- 
ever, the results shoved that a gap between model tips and the end plates 
or increasing the end-plate size for a solid model could result in a 
spirally stable configuration. 


INTRODUCTION 


At the request of the U. S. Aniiy Chemical Corps, a low-speed wind- 
tunnel investigation of the Cook Research Laboratories vortex gliders 
(designated as the Array Chemical Corps model E112 bomblets) was made in 
the Langley stability tunnel. The vortex glider or bomblet was designed 
for use in certain military applications where dispersal of small units 
over very large areas is desirable. The vortex glider under considera- 
tion was designed to rotate rapidly about its spanwise axis while falling. 
This rotation induces a circulation about the glider; consequently, a 
lift is produced. Lift-drag ratios from 1.0 to 2.0 are practicably achiev- 
able and are high enough to be of military value. The advantage of this 
glider over conventional wing-alone or wing-body gliders arises as a 
result of the rapid lift-producing rotation which introduces a gyroscopic 
stiffness that assists in preserving the initial heading. Drop tests by 
Cook Reseairch Laboratories of several configurations of the vortex glider, 
however, showed that the gliders were spirally unstable. 
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The tests of this investigation were made to aid in the development 
of the vortex glider, particularly, to determine a spirally stable con- 
figuration. Static and dynamic lateral stability derivatives were deter- 
mined for the original models and for the models as modified during the 
course of the investigation. In order to expedite distribution of the 
results of this investigation, a detailed analysis has not been made. 


SYMBOLS 


The data are presented about the stability axes and positive direc- 
tion of forces and moments are shown in figure 1. The coefficients and 
symbols used are defined as follows: 

lift coefficient, L/nS 
drag coefficient, D/qS 


cy 

Cl 

Cn 

L 

D 

F 

Y 

% 

M7 


■-x 


lateral-force coefficient, Py/qS 
rolling-moment coefficient, M^/qSb 
yawing-moment coefficient, M 2 ,/qSb 


lift 


drag 

lateral force 
rolling moment 
yawing moment 
dynamic pressure 


moment of inertia coefficient about X-axis, 



Iy 


moment of inertia coefficient about Y-axis, 
radius of gyration about X-axis 
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ky radius of gyration about Y-axis 


^ = 

m 

P 
b 
c 
R 


in 

pSb 

mass 

mass density of air 

span 

chord 


11 = 

n 

V 


2rtnc 

V 

rotational speed of model, revolutions per sec 
free-stream velocity 
yawing angular velocity 
rolling angular velocity 


i*Td 

— yawing-velocity parameter 
2V 


pb 

— rolling- velocity parameter 
2V 

3 angle of sideslip 

7 glide path angle, positive in climb 

dCi 

% " ^ 


•• • ••• 


Sc 


n 


'^R ~ ^ 

p Sp 


c 

2V 

Sc, 

C, = — - 
2V 


Sn„ 

Cn_ = 

2V 

C - ^ 
"^P ~ pb 

S^ 

2V 


SCi 

2V 




S^ 

2V 
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APPARATUS AND MODELS 


All tests were made In the Langley stability tunnel. The static 
yawing and stability derivatives were obtained in the 6- by 6-foot curved- 
flow test section in which curved flight is simulated by causing air to 
flow in a curved path about a fixed model. The rolling derivatives were 
obtained in the 6- foot -diameter test section in which rolling flight is 
simulated by twisting the air about- a fixed model. 

Three models were supplied by Cook Research Laboratories (see fig. 2). 
Models 1 and 2 were mainly made of mahogany with metal end plates. Model 5 
was made of magnesium and provision was made for changing the end plates. 
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Removable tip blocks were supplied with model 3 so that tests could be 
made with different gaps between end plates and model tips. The geo- 
metric characteristics of the models are given in figure 3- The dimen- 
sions of different end plates tested with model 3 are given in figure ^4-. 

A 3A-inch shaft mounted in ball bearings passed throi:igh each model and 
was fastened to a yoke support (fig. 2(b)) which in t\urn was mounted on 
a six-component mecheinical balance system. E>uring the tests the models 
were free to rotate about the shaft and the rotational speed was measured 
with a Strobotac. As a matter of interest it should be pointed out that 
the tunnel models were about three times the size of a practical flight 
model. 


CORRECTIONS 


No jet boundary corrections have been applied to the data. Support 
strut tare corrections, however, have been applied and the lateral-force 
coefficient has been corrected for the buoyancy effect of the static- 
pressure gradient associated with cvu’ved flow. 


TESTS 


Tests of models 1 and 2 were made at a dynamic pressure of 8 and 
l6 pounds per square foot. All tests of model 3 were made at a dynamic 
pressure of 8 pounds per square foot. The airspeeds and Reynolds number 
(based on chord) for the models are given in the following table. 


Model 

q 

V, ft/sec 

Reynolds number 

1 

8 

84 

400,000 

1 

l6 

117 

555,000 

2 

8 

84 

620,000 

2 

l6 

117 

866,000 

5 

8 

84 

436,000 


PRESENTATION OF RESULTS 


The results of the investigation are presented in figures 5 to 48. 
All data presented for models 1, 2, and 3 are based on the following geo- 
metric characteristics. 



Model 

Chord, ft 

Span, ft 

Area, sq ft 

1 

0.750 

2.500 

1.875 

2 

1.167 

2.500 

2.916 

3 

.817 

2.492 

2.036 


The values of Cn.p> presented in table 1 are the 

slopes of the coefficients against 3 obtained near p = 0. 

The derivatives Cn^, Cnp, and C^p presented 

for p other tnan zero were obtained from plots of the coefficients 
against rb/2V and pb/2V similar to those presented for p = 0. Table 2 
is provided as an index to the figiares. 


DISCUSSION 


In order to determine which of the stability derivatives are impor- 
tant for dynamic stability of the vortex glider a study of the coeffi- 
cients of the roots of the stability quartic for the vortex glider is 
necessary. It is well known that the necessary and sufficient conditions 
for stability are that the real parts of the roots of the stability quartic 
be negative. In order that the real parts of the roots shall all be neg- 
ative, it is necessary that B, C, D, E and (BCD-D^-B^D) each be posi- 
tive. The equations for these coefficients have been determined for the 
vortex glider b^ Mr. R. A. Fredette of Cook Laboratories in an unpublished 
report titled, ^'Dynamic Stability of Vortex Gliders." Since the original 
models were spirally unstable and since the algebraic sign of the coef- 
ficient E determines spiral stability or instability, the coefficient E 
was examined to determine which of its terms was of most importance. The 
study showed that the expression for E is similar to that for normal air- 


plane configurations and that for the vortex glider the derivatives 



and Cnp 


are predominant in determining whether E 


is positive or nega- 


tive. The study also showed that the dynamic stability derivatives which 
are found in the expression for E coiild be neglected without altering 
the conclusions with regard to spiral stability.. This resulted in the 
following approximate expression for E 









where k is composed of spin ratio, lift, mass and inertia terms. 
Equating E = 0 to determine the spiral stability boundary eliminates k 
and results in the following expression 


y = 



and 


— — = tan y - - — 


The slope of the spiral stability boundary is, therefore, the negative 
reciprocal of the lift-drag ratio. The approximate E = 0 boundary 
provided a simple criterion for determining whether a configuration was 
spirally stable or unstable and whether the modifications tried were 
beneficial or not. 


The values of Cr 


obtained for model 1 indicate a spirally 
The values of Cr 


-np 
boundary. 


.np and 

unstable configuration as can be seen in figure 45- 
and Cip fell in the unstable region as defined by the E = 0 

This configuration was one of the configurations tested in flight and 
found to be spirally unstable. Although the values of Cn^ and 

obtained for model 2 indicate a spirally stable model at q = l6 as can 
be seen in figure 46, the values for q = 8 indicate spiral instability 
althovigh this is not shown on the figinre. Most of the configurations 
tested were found to have values of Cnp and such that an unstable 

configuration resulted. However, it was found that a gap between the 
model tip and the end plate would result in stability for some of the 
config\arations tested. This is shown in figure 47* The D/L ratio for 
the config\irations indicated on the figiure varies only from 0.91 to 1.05 
so that for simplicity an average boundary was drawn. The configurations 
shown fall well into the stable region. These configurations however have 
a lower l/d ratio than the solid models. This ratio is between 0-95 
and 1.20 for models with a gap. See table 1. 


A spirally stable configuration also resulted for a model with very 
large end plates (end plates number 6, see fig. 4(b)). The bovmdaries 
for this model both with and without a gap are shown on figure 48 and 
both configurations are seen to fall into the stable region. The 
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L/D ratio for the solid model is seen from table 1 to be 2.17 while the 

T 

l/d ratio for the model with 2— — inch gap is 0.95- 

l6 

The inertia of the model with large end plates is probably high and 
the inertia of the actual flight model would also probably be relatively 
high and would require a longer time before full rotation and final glide 
path angle is reached. 


CONCLUDING REMARKS 


The results of the investigation show that a gap between model tips 
and end plates makes some configurations spirally stable, although the 
lift-drag ratio is decreased from that for a solid model. Large area 
end plates ( 19 -inch diameter) result in spiral stability even for a solid 
model with a lift-drag ratio of about 2.0. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 15, 1955- 



William Letko 

Aeronautical Research Scientist 


Approved : 





James L. Williams 
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Thomas A. Harris 

Chief of Stability Research Division 
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TABLE 2.- INDEX OF FIGURES 


Figure 

Model 

Gap size 

End plate 

Driving vane 

1 

Data presented 

5 

1 

0 

See figure 5 

0 

n 

8,i6 

Ci,, Cj), n against p 

6 

1 

0 

See figure 3 



8,16 

Cn< Cl, Cy against p 

7 

1 

0 

See figure 3 



8,l6 

pb 

Cn< Cj, Cy, n against — 

8 

1 

0 

See figure 3 



8,l6 

Cnp. Cip, against p 

9 

2 

0 

See figure 3 



8,l6 

CL) Cp, n against p 

10 

2 

0 

See figure 5 



6, l6 

Cn^ Cy against p 

11 

2 

0 

See figure 3 



8,l6 

rb 

Cn< Cj, Cy, n against — 

12 

2 

0 

See figure 3 



8,l6 

Cjj, Cj, Cy, n against ^ 

13 

2 

0 

See figure 5 



8,l6 

Cnj.< P 

14 

2 

0 

See figure 3 



8,l6 

p 

15 

3 

0 S ll 2-1 

8’ 32’ 16 

1 



8 

Cp, Cp, n against p 

16 

3 

0 ^ ll 2l 
8 52 16 

1 



8 

Cn; Cj, Cy against p 

17 

3 

7 

0, 2-5- 
l6 

1 



8 

rh 

Cn^ Ct, Cys D against • — 
^ 2V 

18 

3 

0, 2I 

16 

1 



8 

Ci^, Cy^ against p 

19 

3 

0, 2^ 
16 

2 



8 

Cp, n against p 

20 

3 

0, 2-2- 
16 

2 



8 

Cq, Cj, Cy against p 

21 

3 

0, 2-L 

16 

* 3 



8 

Cl< Cp, n against _ p 

22 

3 

0, 2-L 
16 

^ 3 



8 

Cn. Cl, Cy against p 

23 

3 

0, 2-^ 
16 

“ 3 



8 

rb 

Cn» Cl, Cy, n against — 

24 

3 

16 

® 3 

> 

r- 

f 

8 

Cj^, Cy^ against p 


®Long dimension perpendiciilar to chord. 
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TABJ.E 2.- INDEX OF FIGURES - Concluded 


Figure Model Gap size End plate Driving vane 


Data presented 


o C£|, n against p 

6 Cn^ C^, Cy against P 

6 Cl, Cl, n against p 

8 Cn, C^, Cy against p 

® Cn, C 2 , Cy, n against 


8 Cn^. Ci^, Cy^ against p 


8 Cj^, Cj), n against p 
8 Cn, Cl, Cy against p 
8 Cjj, Cj, Cy, n against 


8 Cnj,, Ci^, Oy^. against p 


8 Cl> Cj), n against p 
8 Cn, C;, Cy against p 
8 Cn, Cl, Cy, n against 


8 Cn^., Cij., Cy^, against p 


Cl> 

Cd, 

n 

against 

p 

Cn, 

Cl, 

Cy 

against 

p 

Cl. 

Cd, 

n 

against 

p 

Cn> 

Cl, 

Cy 

against 

p 

Cl' 

Cj), 

n 

against 

p 

Cn> 

Cl, 

Cy 

against 

p 


figure 3 


Spiral statility boundary 
6 Spiral stability boundary 
8 Spiral stability boundary 

8 Spiral stability boundary 


®-Long dimension perpendicular to chord. 
'^Long dimension parallel to chord. 


^1 OJ 
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(a) Models 1 and 2. 


L- 8871^.0 



Figure 2.- Photographs of models. 
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(b) Model 5- L- 9076 I 4 . 


Figure 2.- Concluded. 
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End plate 4 


End plate 5 


End plates 4,5,6 thick 

(b) End plates 4, 5 , and 6, 
Figure Continued. 
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End view 
End plate 7 

( c ) End plate 7 • 
Figure 4.- Concluded. 
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Figure 6.- Variation of Cj, and Cy with angle of sideslip p. 

Model 1; q = 8 and l6. 
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-20 -16 -12 -8 -4 0 4 8 12 16 20 

Angle of sideslip, J9, deg 


Cjn } C 7 1 ajid. Cv 

T 

Model 1; q = 8 and l6. 


Figure 8.- Variation of 


with angle of sideslip 3 



NACA RM SL55J26 


• » ••• •« 


mmm 

PPPPPPPP 


M 

PPbPBI 

m 

pp 

PI 



jSass 

safaassjgfe^cassigs 

S 

sa 

s^sasa^ual 

tss 

saps 

san 

aass 

Pi 

Bass 

IPPPPPP 

SBasasaassasa 

P 

sa 

PPPPPPi 

caaaaaa3a&as:l 


Q 

O 8 
□ 16 












g^s~!!«SSae!ssSaSS^‘®^gaa»SKn«25ii«S 
astatssgsaBSists tasssssssi 


rps 


__ ~ »w»ns««t»5w»§»^»{?»»S»^»»*snia 

ia^ ^a ^^@aga ^^pp giggs gasaP ag BPBP pppsi!uaigss^S5^gg^i 

asg^S^fflissa^aa^Bgsagasa^Ba^^iajnaSS^Bsgniggginjnn jig ggjaggiSssj^^ 


SB B I SU iilflnmij SSmSmS SSIZBSS AB^tBB B iilSISi lB S B UxB BltiSSSBiBBi mm A i S M SS Uiilt?«BiBBB»*aSmiB*«BSB£SKS mMlSlt* 

^ gs:aa|pS ji^^^j^ |a«5S^{gaaffl as as; g;: 


gwatsatstrisaataKsasst! 


^1:S!Ju!aS?Sa_ — 

gBBBBBaSsBHjgjga g 


|gsgagai|Bgs|s|§|is^ jsjlilssy^SpB 


^ss;SS53 ^sss Bst^ ^ssass 
SSaaJisa ssi^»^5sa;^s ^EataSH 






SI Silaata aSgasa 


: ga; t g;;»s ; sa;;a:!aa ga8»»ta;;aa;;! 
Si~»~Spg|n^SiSi 


B gBjjjjjBgHjjHSBjHIgHBljaSjmisnB s n g 


iPPPPPPPPPPlPPPr 


5B5g!^^BgaBBBBBBsaiBPPPPBB Bte 5-!ggBBBigBaBiaBBggiagssiBI 

pBPPPPPPPBPPPPPpPPppfiP?|:PPp||]||PP|i|ipPP||||ipi; 

^§p§p^i^p^pBP^sag aP P P PPBlalSPPPBSi.Sg:PnaPPPPPPP~ii 

gm^m^^iiIIjnni!^iaa^jBiSi‘i^P^SSgag:5Sgg a^^g BByjgiHgt8gggiygiiyggigigj g;;gHSK!agSgSIS5 HiaII! ll! II 




a jgjigfesBPBPBBBBB BB BSgagagiss Bgr 

ppppppipippipipppppppp 


[■iBigilBBBPBiss 



PbPPPPPPPsPPP BBPPPPBPP 




ppppipppsps^rstiiatil^rsippp 

IPPPPtiPPSP 


sSgaS PPPssPBPPgg 

g pB PPmgBPBsgPP^ 


isiSgaiE-igggHiglllBBBc 

ipyg^i|BgBlp!!!8glJll! 


ii!gi3n:iHa=- - 


SBpPgglPg 

intisnsntiiisa s 
ga^|a:aata: .. 


Cr 


IPii: 
iP®JPP“ 


IPPigli!! ~ ’■ II 

^s«iBs;a .. F 

:pEtu::nssm!^» 

SggSpiili ' ' 




Hi! fiPPPPlfJl 

Iniin: pBl^l^^sjijijj!!'^ 


n^SHSSfSsaaSiSHr 
a :s;; ;ts; ;;:n:ls; Ss;;u^s 




gBB 


gppiaigii 


1 

ppI 


PPI 


6 20 


Angle of sideslip, /3, deg 


Figiire 9*- Variation of Cl, Cl, and n with angle of sideslip p. 

Model 2; q = 8 and l6. 









Figure 11.- Variation of C^f an^ n with rh/2V for 3 = 0. 

Model 2; q = 8 and l6. 
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Figure I3.- Variation of . C^, , and Cy with angle of sideslip p. 

Model 2; q = 8 and 16. 



-20 -16 -12 -8 -4 0 4 8 12 16 20 

Angle of sideslip, j3, deg 

Figure l4.- Variation of > C 7 , and Cv with angle of sideslip 3. 

T p p 

Model 2; q = 8 and 16 . 
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Angle of sideslip, /9, deg 


Figure 15.- Variation of C-j^, Cp, and n with angle of sideslip p 

for model 3 with different size gaps. End plate 1> driving vane on; 

= 8 . 
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Angle of sideslip, j8, deg 


Figure l6.- Variation of Cn^ and Cy with angle of sideslip P 

for model 5 with different size gaps. End plate 1; driving vane on; 

q = 8. 
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Figure 18.- Variation of Cj, ^ C 7 , and Cy with angle of sideslip p 

r "^r 

for model 5 with 0- and 2^ - inch gap size. End plate 1) driving vane 
on; q = 8. 
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Figure I 9 .- Variation of C-^, C-q, and n with angle of sideslip P 
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Figure 20.- Variation of Cj^, Cj, and Cy vith angle of sideslip p 

for model 5 vith 0- and 2-i - inch gap size. End plate 2; driving 

16 

vane on: q = 8. 
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Angle of sideslip, /3, deg 


Figure 21.- Variation of C^,, and n with angle of sideslip p 

7 

for model 5 with 0- and 2-^ - inch gap size. End plate 3 with long 
dimension perpendicular to chord; driving vane on; q = 8. 
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Angle of sideslip, fi, deg 


Figure 22.- Variation of Cq, and with angle of sideslip p 

7 

for model 3 with 0- and 2-1— - inch gap size. End plate 3 with long 

16 

dimension perpendicular to chord; driving vane on; q = 8. 
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Figure 25.- Variation of Cn, Cy^ and n with rh/2V for model 5 

7 

with 0- and 2.-^ - inch gap size. End plate 5 with long dimension 
Id 

perpendicular to chordj driving vane on; q = 8. 
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Fig\are 2k.- Variation of and Cy^ with angle of sideslip 3 

7 

for model 5 with 0- and 2-!—- inch gap size. End plate 3 with long 

16 

dimension perpendicular to chord; driving vane on; q = 8. 
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Figure 25-- Variation of Cj^, Cp, and n with angle of sideslip p 

for model 3 with 0- and '2^- - inch gap size. End plate 3 with long 

Id 

dimension parallel to chord; driving vane on; q = 8. 
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Figure 26.- V^iation of Cn, Cj, and Cy with angle of sideslip p 

7 

for model 3 with 0- and 2 — - inch gap size. End plate 3 with long 

16 

dimension parallel to chord; driving vane on; q = 8. 
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Figure 27 .- Variation of and n with angle of sideslip p 

7 

for model 3 with 0- and 2— inch gap size. End plate 4; driving 

16 

vane on; = 8 . 



NACA RM SL55J26 




Cy O 


08% 

I 

.04% 
O ^ 
-.04 
-.08 
-.12 
■08 
04 

r> 

O 

-.04 

-08 


mmMmmmmmmmmmmi 


Gap size 
O O 
° 



BBSqpa 


sr;<a8g8sa;f^ 



BB 

BBBgi 


SSi 

S@Bfifeg=,Sr" 






iHoniBss:;; 




5 55^caS n ? n3 CT ^i^}|i^ i^ ^ Eii!ini 5:::S 


IBSSBBal 


BSBB§ 




«« p«225 ^^9 £SS3 S5S 

8 ^S5 


gii.-.-- -| 






"IBHllBSBBSBI 


IBaBSBSSs 


IBBBaaBpBEsB 


aWBBWBBBBBBBBB BBBBBBBBBBBIgg i 


g ^ ^ :»pg; 8s; s» 55 






sH^garpl^jp l^i 






-20 -lb -!2 -tf -4 0 4 8 12 16 20 

Angle of sideslip, j3, deg 

Figure 28.- Variation of C 7 , and Cp with angle of sideslip 3 

Y 

for model 5 with 0- and 2-~ -inch gap size. End plate driving 
vane on; q = 8 . 
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Figure 30.- Variation of 

for model 3 with 0- and 2— 

16 

on; q = 8. 


ajid Cy^ with angle of sideslip p 
- inch gap size. End plate 4; driving vane 
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Figure 32.- Variation of Cj, and Gy with angle of sideslip 3 

for model 5 with 0- and 2^^ - inch gap size. End plate 5; driving 
vane on; q = 8. 





Figure 33.- Variation of Cn^ C^, Cy, and n with rh/2V for model 3 
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Figure 5^-- Variation of Cn^, and Cy^ with angle of sideslip 3 


for model 5 with 0- and 2-^ - inch gap size. End plate driving vane 

Id 

on; q = 8. 
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Figure 35.- Variation of Cl, Cj), and n with angle of sideslip p 

for model 3 with 0- and 2^ _ inch gap size. End plate 6; driving 
vane on; q = 8. 
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Figure 36.- Variation of and Cy with angle of sideslip p 

for model 3 with 0 - and 2 ^ -inch gap size. End plate 6; driving 


vane on: q = 8. 
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Figure 38.- Variation of 0^^,, and Cy^, with angle of sideslip p 

7 

for model 3 with 0- and 2-1— - inch gap size. End plate 6j driving vane 

16 

on; q. = 8. 
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Figure 59-- Variation of Cjj, and n with angle of sideslip p for 

7 

model 3 with 0- and - inch gap size. End plate 7j driving vane on; 
q = 







model 3 with 0- and 2^ - inch gap size. End plate 7; driving vane on;- 
q = 8. 
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j Figure 4l.- Variation of Cl, Cl, and n with angle of sideslip P for 

i model 3 with 2-i- inch gap size. End plate 1; driving vane off; q = 8 
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Figure 42.- Variation of and Cy with angle of sideslip p for 

Y 

model 3 with 2-i- - inch gap size. End plate 1; driving vane off; q = 8. 
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Figure 43.- Variation of Cj^, Cn, and n with angle of sideslip 3 for 

model 3 with end plates off and with end plates and 2-^ - inch filler 

l6 

blocks off. Driving vane on: q = 8. 
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Figure 44.- Variation of 0^, Ci, and Cy with angle of sideslip 3 for 
model 5 with end plates off and with end plates and 2-2- - inch filler 

o 16 

blocks off. Driving vane on: q = 8. 
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Figtire 45.- Spiral, stability boiindary for model 1 with position of model 1 

with respect to the boundary, q = 8. 




NACA m SL 55 J 26 



.007 

.006 

.005 

.004 

.003 

.002 

. 00 ! 

% 0 
-. 00 ! 
-.002 
-.003 
-.004 
-.005 
-.006 



-.007 -.006 -.005 -.004 -.003 -.002 -. 00 ! 0 00 ! 002 003 



Figure 46.- Spiral stability bovmdary for model 2 with position of model 2 

with respect to the boundary, q = l6. 
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Figure 47.- Spiral stability boundary for several models and position of 
models with respect to the boundary, q = 8. 




